How to detect level crossings without looking at the spectrum 
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We remind the reader that it is possible to tell if two or more eigenvalues of a matrix are equal, 
without calculating the eigenvalues. We then use this property to detect (avoided) crossings in the 
spectra of quantum Hamiltonians representable by matrices. This approach provides a pedagogical 
introduction to (avoided) crossings, is capable of handling realistic Hamiltonians analytically, and 
offers a way to visualize crossings which is sometimes superior to that provided by the spectrum. We 
illustrate the method using the Breit-Rabi Hamiltonian to describe the hyperfine-Zeeman structure 
of the ground state hydrogen atom in a uniform magnetic field. 
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I. INTRODUCTION 

Crossings and avoided crossings occur in the spectra 
of many physical systems such as atoms ji molecules j 2 - 
semiconductors^ and microwave cavities^ They can oc- 
cur as a result of tuning an external (e.g. electric) fielo-i 
or as the consequence of varying an internal (e.g. in- 
ternuclear) coordinate^ 6 - for instance. A number of in- 
teresting physical phenomena have been associated with 
such (avoided) crossings. For example an eigenstate 
transported in a closed circuit around a crossing in the 
spectrum picks up a Berry phased As another example, 
points in the spectrum where a large number of eigen- 
values cross correspond to 'hidden' symmetries of the 
physical system.— These symmetries are hidden in the 
sense that they are not evident a priori as observables 
that commute with the Hamiltonian. As a third example 
eigenvalue avoidance in the spectrum signals the emer- 
gence of quantum chaos in a physical system^ 

Inspite of its importance as a basic phenomenon ubiq- 
uitous in physics, few introductory texts treat (avoided) 
crossing in any detail. Some questions that arise in the 
context of simple physical systems, and that could be 
considered at the (under)graduate level, are : 

1. Is there a way to predict the total number of 
(avoided) crossings in the spectrum ? 

2. Is there a systematic way to locate all the (avoided) 
crossings in the spectrum ? 

3. Is there a way to identify the physical mechanisms 
responsible for the occurrence of (avoided) cross- 
ings in the spectrum ? 

4. Can the degeneracies - if any - in the spectrum be 
thought of as crossings ? 

5. Sometimes the crossings arc hard to discern in the 
spectrum - can they be visualized in a clearer way 

? 



In this article we employ an algebraic method that ad- 
dresses these basic questions as well as some others. It 
allows for a simple but systematic approach to (avoided) 
crossings. In describing this approach we reintroduce to 
the study of (avoided) crossings a very useful but perhaps 
neglected mathematical tool, the discriminant i 10 : 11 : 12 : 13 
In a series of articles ) 14 ' 15 : 16 we have demonstrated in de- 
tail how use of algebraic techniques is a powerful way 
to locate (avoided) crossings in the spectra of quantum 
mechanical systems. We have shown that they are not 
only capable of locating (avoided) crossings without re- 
quiring solution of the Hamiltonian - a fact well known 
to mathematicians and not unknown to physicists - they 
can also find (avoided) crossings when the Hamiltonian is 
not completely determined^ As one example, algebraic 
techniques allow us to derive a new class of invariants 
of the Breit-Rabi Hamiltonian. These invariants encode 
complete information about the parametric symmetries 
of the Hamiltonian^ As another example, the use of al- 
gebraic methods allows us to detect the breakdown of the 
Born-Oppcnhcimer approximation for molecules assum- 
ing only that the complicated molecular potentials are 
realJ£ 

Although in the work just mentioned algebraic tech- 
niques were used in the context of advanced research, 
namely Fcshbach resonances, we show below that they 
also form an effective pedagogical tool. Indeed we believe 
the exposition in this article could easily be included in 
the physics curriculum for (under)graduates strik- 
ing demonstration that enhances their understanding of 
quantum mechanics as well as linear algebra. 

The rest of the paper is arranged as follows. Section 
II contains a simple mathematical introduction using a 
2x2 matrix. Section III generalizes this to the case of an 
n x n matrix . Section IV demonstrates the technique 
developed so far on the ground state hydrogen atom in 
a uniform magnetic field. Section V suggests some exer- 
cises for the reader; Section VI supplies a discussion. 



II. A 2 x 2 MATRIX 



III. AN n x n MATRIX 



In order to motivate the general case we first consider 
a real-symmetric 2x2 matrix 



H(P) 



Ei V 
V E 2 



(1) 



with (unknown) eigenvalues Ai. 2 . The notation implies 
that all the matrix elements depend on some tunable pa- 
rameter P. Also E\ : 2 could be the bare energies of a 
two-level quantum system, which are mixed by the per- 
turbation V. Usually to find Ai j2 we solve the equation 

\H(P) - A| = 0, (2) 

where A is a parameter. Eqs.fl]) and © yield 

A 2 + dA + C = (3) 

where Co = E\Ei — V 2 and C\ = —(E\ + E 2 ). However 
Ai ; 2 also satisfy Eq.©, i.e. 

(A-A!)(A-A 2 ) = 0. (4) 

Comparing the coefficients of Eqs.Q and we find 



Co — A1A2, 

Ci = -(Ai + A 2 ). 



(5) 



We are interested in crossings in the spectrum of H(P) 
and therefore consider the discriminant^ defined by 



A = (Ai-A 2 ) 



(6) 



A little tinkering with Eqs.© shows that the discrimi- 
nant can be re-written solely in terms of the coefficients 
of Eq.© : 



A = C\ -4C . 



(7) 



It is important to note that we did not actually calculate 
A12 at any point in the discussion so far. Clearly, A = 
whenever a level crossing occurs in the spectrum of H(P). 
For example, if E x = E 2 = P, and V = 2P, then A = 
4P 2 , and there is a level crossing at P = 0. This may 
be verified by calculating the eigenvalues Ai j2 = -P, 3P. 
Note that a single crossing in the spectrum corresponds 
to a double root of the discriminant. 

We see from this example that use of the discriminant 
transforms the problem of finding crossings in the spec- 
trum to a polynomial root-finding problem. Further, it 
enables us to avoid calculating the eigenvalues. Lastly, 
it provides the locations of all the crossings in the spec- 
trum. In the next section we generalize these statements 
to the case of an n x n matrix. 



For a real symmetric n x n matrix H{P) all of whose 
entries are polynomials in P the characteristic polyno- 
mial is 



\H{P) - A| = V C m A 



m=0 



(8) 



where the coefficients C m are all real. The discriminant 
is defined &s~ 



DlHiP)}^]]^-^) 2 , 



(9) 



i<j 



in terms of the n eigenvalues A^. It can also be written 
purely in terms of the n + 1 coefficients C m >i£ 



D[H(P)] = 



Cn Cn — 1 



Co 



Co 



nC„ 



nC n 



Ci 



2C 2 Ci 

(10) 

In practice we will calculate discriminants of characteris- 
tic polynomials using the built-in Discriminant func- 
tion in Mathematical or Maple. In addition we will use 
the following 'toolbox' of general results in our discussion 
of (avoided) crossings below : 

1. The real roots of D[H(P)] correspond, as in Section 
Hn to crossings in the spectrum of H{P). 

2. The real parts of complex roots of D[H(P)] cor- 
respond to avoided crossings in the spectrum of 
H(P). For a proof, see RefFlfJ. 

3. A crossing is defined as the intersection of two 
eigenvalues. Hence the simultaneous intersection 

m 



of m eigenvalues gives rise to 
crossings. 



= m(m-l)/2 



4. Every (avoided) crossing contributes, as in Section 
ITu a factor quadratic in P to D[H(P)]. For a full 
proof, see RefflEI. 

5. Since H(P) is real symmetric, the eigenvalues Aj 
are all real. It follows from Eq.© that D[H(P)\ > 
0. Hence Log[D[H (P)] + 1] > and goes to zero at 
every crossing. We will plot this function in order 
to visualize crossings. 
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IV. THE HYDROGEN ATOM 

A simple but real example of a physical system whose 
spectrum exhibits both crossings and avoided crossings 
is a ground state hydrogen atom in a uniform magnetic 
fields Such an atom is accurately described by the Breit- 
Rabi HamiltonianSI 



We now systematically locate all the (avoided) cross- 
ings in the spectrum. There is a 6-fold real root at B = 0, 
which points to the only crossings in the spectrum (Scc- 
tion lllll l.). Conventionally these 'zero-field' crossings are 
called degeneracies and are not considered to be cross- 
ings. Hence we will maintain that there are no crossings 
in the spectrum of H]^ . Eq. (Q3J) has two complex roots: 



H BR = AL-S + B(aS z + bI z ), 



(11) 



where / and S indicate the proton and electron spin 
respectively, and B is the magnetic field along the z- 
axis. A equals the hyperfine splitting and a = g e ^B and 
b = ,9 p /ijv where g e ( P ) are the electron(proton) gyromag- 
netic ratios and Hb(n) are the Bohr (nuclear) magnetons 
respectively. The numerical values for these constants 
were obtained from RefH^ and are to be found in the fig- 
ure captions in this article. In the basis |Mj, Ms) which 
denotes the projections of / and S along B, the states 
are 



I- -) I- -) I- 

'2' 2 M 2 ' 2 M 2' 



-1 



>,| 



1 -1 



(12) 



In this basis, the representation of Eq. (fTT|) is 



-A + 2(a - 
2A 



■ b)B 



In order to calculate discriminants the basis (|12|) is all 
we need but we would also like to correlate our results 
to the spectrum, where we will use the basis \F,Mp) 
to label the eigenstates of (fT3|) . Here F — I + S is the 
total angular momentum and Mp its projection along the 
magnetic field. This labelling is loose as strictly speaking 
in the presence of a magnetic field M p is a good quantum 
number but F is not. 



B 



a 



(15) 



This conjugate pair corresponds to a single avoided cross- 
ing at B = (Sectioning, and 4.). From Ea.((I5|l we can 
see that A, the hyperfine coupling of the two spins, is the 
physical mechanism responsible for the avoided crossing. 
Setting A = turns the avoided crossing into a zero-field 
degeneracy. 

We can confirm the predictions of the discriminant 
(Ecim]) by looking at FigOJa). At B = 0, three states 

in the upper manifold coincide, giving rise to ( ^ 

crossings (Section |Hl] 3.). Further, each crossing con- 
tributes a factor quadratic in B (Section IIIII 4.) to the 
discriminant (Eq fT4|) , which therefore should - and does - 
contain a factor of B 6 . This shows that the discriminant 
accounts for the (hyperfine) degeneracies of the spectrum 
by counting them as crossings. 

An avoided crossing occurs at B = between the states 
drawn with dotted lines. In Fig[ljb) a logarithmic rep- 
resentation of D[H b jf^] shows a dip corresponding to the 
degeneracies at B = 0. 



The 6/0 case 



When b ^ 0, we find from Eq.JTS 



A. The b = case 

To begin we recollect that \b\ \a\, since the proton 
is more massive than the electron.— We set b = in 
Eq. tfTB"]) . calculate its characteristic polynomial, and find 
the discriminant to be 



1 



D[H B =°] = Y^A 4 a 6 B 6 (A 2 + a 2 B 2 ) . 



(14) 



Considering D[H b jf^] as a polynomial in B, the total 
number of its roots equals its degree, namely 8. Using 
the results 1., 2. and 4. from Section ITTT1 we can then 
say that the total number of the avoided crossings plus 
the number of crossings in the spectrum of is ex- 

actly 8/2=4, a fact verified below. Thus the discriminant 
allows us to predict the number of (avoided) crossings be- 
fore they are actually found. 



D[H BR ]= ±(a + b) 2 B e [A 2 + (a-b) 2 B 2 ] 2 

x [(a + b)A + 2abB] 2 [(a + b)A - 2abB] 



(16) 



The 6-fold real root at B = persists from the b = 
case (EqQjJ, but now there are ('authentic') crossings 
for B ^ at 



B = ± 



(o + b)A 
2ab ' 



(17) 



'Switching' b off and on therefore reveals the physical 
mechanism behind the appearance of crossings in the 
Breit-Rabi spectrum : it is the interaction of the pro- 
ton spin with the external magnetic field B. 
The complex roots of Eq. (fl"6]) are 



B 



a — o 



(18) 
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FIG. 1: (a) The spectrum of H B ~ R a for atomic hydrogen (S = 
1/2,1 = 1/2) with states labelled using the \F,Mf) basis. 
The parameters used were A — 0. 0473cm -1 and a — 9.35 x 
10 _8 cm _ G _1 . The spectrum was calculated by analytically 
diagonalizing the matrix in Eq. (|13p and numerically plotting 
the eigenvalues, (b) liO%[N a D[H b B = R ] + 1] plotted using the 
same parameters as in (a) and the scaling factor No = 10 3 ° 
used to optimize visibility. The overall logarithmic trend is 
clearly visible. 
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and imply an avoided crossing at B = as in Section 

EZSK EqgB). 

We can confirm the predictions of the discriminant (Til?]) 
by looking at Fig{2j The crossings which occur at ±16.6 
T (from Eq |17[) for actual values of the hydrogenic param- 
eters a, b, and A cannot be seen in the spectrum [Figl^a)] 
as the energy level pairs are separated by less than the 
width of the lines used to plot them, a p oint made ear- 
lier in this journal (see Fig. 2 in Refr20l). However the 
logarithmic representation [FigJ2{b)] of the discriminant 
on the same scale clearly shows dips at all the crossings. 
A scaled spectrum has been shown in Figf2^c) using a 
much larger relative value of b in order to display the 
crossings clearly. This example illustrates how the dis- 
criminant can sometimes prove superior to the spectrum 
in displaying crossings. For another example see ReffiH. 



FIG. 2: (a) The spectrum of Hbr for atomic hydrogen (S = 
1/2,1 = 1/2). The parameters used were A = 0.0473cm" 1 , 
a = 9.35 x 10" s cm" 1 G" 1 and b = -1.4202 x 10~ 7 cm _1 G -1 . 
The spectrum was calculated by analytically diagonalizing the 
matrix in Eq. (|13[) and numerically plotting the eigenvalues. 
The crossings at B — ±16.6 T cannot be resolved on this scale, 
(b) Log[N b D[HBR] + 1] plotted using the same parameters as 
in (a) and the scaling factor Nb = 10 10 . All three crossings 
are distinctly indicated by the minima of the discriminant; 
those corresponding to crossings at B — ±16T actually do 
reach zero, but appear shorter here due to the resolution of 
the graphics. The overall logarithmic trend of the plot is 
evident, (c) The same spectrum as in (a) but redrawn using 
the parameters A = 0.1, a = 0.01, and b — —0.1. The larger 
relative value of b ensures that the crossings at B / can be 
seen in the plot. 
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V. SUGGESTED EXERCISES 

1. Prove that the product in Eq.® contains n(n — 
l)/2 factors. 

2. Justify the presence of A A in Eq. (fT4"|) . 

3. For the parametric symmetry a = b in Eq. (|16[) show 
that there are crossings but no avoided crossings in 
the spectrum of Hbr - 

4. Plot D[Hbr] as a function of B and try to iden- 
tify the zeros that correspond to crossings. This 
is an exercise designed to show that the highly 
nonlinear nature of the discriminant implies that 
each of its terms dominates in a different regime 
of B. This makes it difficult to include all the fea- 
tures of the discriminant on a single scale unless 
a smoother representation - such as a logarithmic 
one - is adopted. 

5. The Wigner von- Neumann non-crossing rule^ says 
"States of the same symmetry (i.e. quantum num- 
ber) do not cross, except accidentally." Verify the 
rule for the \Mp) states in FigsQJa) and^c). That 
is, show that states which (avoid) cross possess 
(same) different Mp's. 

VI. DISCUSSION 

The examples presented above illustrate that the dis- 
criminant is an elegant but simple device for locating and 
counting (avoided) crossings. Further, it is an effective 
tool for investigating the physical mechanisms behind the 
occurrence of (avoided) crossings. Lastly, visualization of 



the discriminant offers an alternative to locating cross- 
ings in the spectrum. It must be noted however that the 
discriminant yields no information about which eigenval- 
ues avoid or intersect, or about the eigenvectors. Also 
shallow avoided crossings do not show up in the logarith- 
mic representation, especially if they are near to cross- 
ings, which give rise to strong features in the discrimi- 
nant. For all such information the spectrum has to be 
consulted. 

The technique we have presented can be used algorith- 
mically on Hamiltonians which are polynomial in some 
parameter P. and which can be represented by finite di- 
mensional matrices. Examples are a spin 1/2 particle 
in a magnetic field^ (the archetypal two-level system) 
and the hydrogen atom in an electric field^ (usually pre- 
sented as an example of degenerate perturbation theory) . 
However the method can also yield insight into physi- 
cal systems whose Hamiltonians are usually truncated 
to a finite dimension for practical calculations such as 
the nucleus modeled as a triaxial rotor^ and a polar 
molecule in an electric fields Another interesting appli- 
cation is the calculation of critical parameters of quan- 
tum systems^ since the critical point occurs at a cross- 
ing. An example using the Yukawa Hamiltonian has been 
treated in Refl27l. A list of physical systems in atomic 
and molecular physics to which algebraic methods can 
be applied is provided in Reffltl 

To conclude we have presented an algebraic technique 
for systematically analysing (avoided) crossings in the 
spectra of physical systems and pointed out its useful- 
ness as a pedagogical device. It is a pleasure to thank 
P. Meystre for support. This work is supported in part 
by the US Army Research Office, NASA, the National 
Science Foundation and the US Office of Naval Research. 
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